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ABSTRACT 


The  large  subsonic  "academic"  wind  tunnel  at  NFS, 
powered  by  two  counter-rotating  fans,  never  achieved  the 
design  specifications.  The  reason  for  the  poor  performance 
of  the  tunnel  was  unknown  but  believed  to  be  due  to  either 
poorly  designed  fan  blades  or  to  separation  in  the  diffuser. 
Being  easier  to  analyse,  the  fan  blades  were  chosen  for 
initial  study. 

The  approach  used  is  a  new  blade  element  method  for 
calculating  the  performance  of  high  and  intermediate 
solidity  fans.  Although  this  new  method  predicts  some  devia¬ 
tions  from  the  original  isolated  blade  analysis,  it  was 
found  that  the  design  was  adequate  and,  therefore,  the 
tunnel  problem  is  most  probably  due  to  separation  in  the 
diffuser. 
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I.  INTRODUCTION 


The  "academic"  wind  tunnel  at  NFS  is  inadequate  for  most 
tasks,  due  to  the  excessive  level  of  turbulence  at  high 
speed.  This  tunnel  was  designed  for  a  potential  speed  of  200 
Knots  in  the  test  section,  but  this  value  has  not  been 
reached  yet. This  project  was  undertaken  with  the  purpose  of 
evaluating,  by  a  different  method  from  that  used  for 
designing,  the  performance  of  the  tunnel  in  order  to 
establish  possible  causes  of  the  problem. 

The  two  counter  -  rotating  fans  could  be  the  cause  of 
the  trouble  since  they  were  designed  using  a  blade  element 
theory,  that  does  not  take  account  of  the  three  dimensional 
effects  that  are  present  at  medium  and  high  solidities. 
These  effects  cause  changes  in  the  apparent  velocity  vector, 
therefore  the  blades  were  analysed  by  a  "  New  Blade  Element 
Method  for  Calculating  the  Performance  of  High  and 
Intermediate  Solidity  Axial  Flow  Fans  "  due  to  Borst 
[Ref.  1:  p.l].  This  method  determines  an  induced  angle  of 
attack  which  changes  the  apparent  velocity  past  the  airfoil 
and  it  has  been  found  to  be  a  reliable  measure  of  the  three 
dimensional  effects. The  data  required  for  this  method  comes 
from  the  two  dimensional  airfoil  data,  that  is  available  in 
the  current  literature.  This  method  was  used  in  reference  1 
(p.lO)  assuming  zero  drag,  and  the  predictions  showed 
excellent  agreement  with  measurements.  The  small  difference 
between  the  results  can  be  attributed  to  the  assumption  of 
zero  drag.  The  performance  of  the  first  and  second  stage 
blades  of  the  academic  wind  tunnel  is  evaluated  assuming  non 
zero  drag,  so  it  is  e.xpected  that  the  results  will  be 
realistic . 

A  possible  cause  of  the  inability  to  reach  the  specified 
speed  could  be  excessive  losses  around  the  tunnel  circuit, 
principally  at  the  first  diffuser  and  the  corners.  The  first 


diffuser  has  a  high  angle  of  divergence  which  can  provide  a 
flow  separation,  further  losses  and  turbulence.  The 
evaluation  of  losses  around  the  tunnel  is  made  by  means  of 
the  method  given  by  RAE  &  POPE  [Ref.  2:  p.87],  which  breaks 
the  tunnel  down  into  sections  and  calculates  the  losses  for 
each  one.  The  four  corner  losses  are  evaluated  by  three 
methods : 

1  considering  one  third  of  the  losses  due  to  skin 
friction  and  two  third  losses  due  to.  rotation;  an 
empirical  relation  given  by  RAE  &  POPE  [Ref.  2:  p.90], 
was  used  to  account  for  the  different  types  of  corner 
vane , 

2  considering  empirical  values  found  by  testing  and 

3  considering  variation  of  resistance  coefficient  with 

tap/ chord  ..ratio  of  thin  corner  vanes  given  by  PANKHURST 
HOLDER  [Ref.  3:  p.93]. 

The  first  diffuser,  since  it  has  variable  divergence, 
was  broken  into  four  parts  and  the  losses  evaluated  part  by 
part  through  the  method  given  by  RAE  &  POPE  [Ref.  2:  p.89], 
this  shows  that  the  losses  are  strongly  dependent  on  the 
diffuser  angle. 

After  all  losses  were  evaluated  ,  a  relation  between  the 
air  inflow  velocity  to  the  blades  and  the  power  output  from 
the  blades  is  found  and  the  non-stall  operation  envelope  is 
determined.  The  temperature  and  pressure  assumed  by  LARSON 
in  the  original  blade  design  [Ref.  4:  p.5],  are  T=  100°  F 
p=  2,246  (Ibf/ft  )  and,  therefore  the  values  of  kinematic 
viscosity  (;/)  of  0.000166  (ft  /  sec)  and  sound  velocity  of 
1160  ft  /sec,  as  presented  in  Appendix  C,  will  be  used  for 
evaluation  of  the  Reynolds  and  Mach  numbers. 


II.  EVALUATION  OF  LOSSES 


A.  INTRODUCTION 

As  mentioned  in  chapter  one, we  will  use  the  procedure 
outlined  by  RAE  &  POPE  [Ref.  2:  p.87],  that  consists  of 

breaking  the  tunnel  down  into  cylindrical  sections, 
expanding  sections,  contracting  sections  and  corners,  in 
order  to  evaluate  the  tunnel  losses. 

The  sections  are  numbered  in  such  way  that  those 
sections  with  the  same  characteristics,  that  is,  those  that 
use  the  same  loss  equation,  are  kept  together.'  This 
facilitates  the  evaluation  by  computer.  Then,  the 
cylindrical  sections  are  numbered  from  1  to  4,  the  expanding 
sections  from  5  to  11,  the  contracting  section  and  spinner 
take  the  number  12  and  13  respectively,  and  finally  the 
corners  from  14  to  17. 

The  Figure  2.1  shows  the  sections  and  their  designated 
numbers,  as  they  are  discriminated  below: 

DISCRIMINATION 

test  section 
fan  duct 

constant  area  duct 
constant  area  duct 
first  diffuser  (sec  1) 

"  diffuser  (sec  2) 

"  diffuser  (sec  3) 

"  diffuser  (sec  4) 

second  diffuser  (sec  1) 

"  diffuser  (sec  2) 

"  diffuser  (sec  3) 
contraction  cone 
spinner 

first  corner  with  full  size  vanes 


SECTION 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


15 

second 

corner 

with 

full 

size 

vanes 

16 

third 

corner 

with 

full 

size 

vanes 

17 

fourth 

corner 

with 

half 

size 

vanes 

The  dimensions  of  the  sections,  which  are  used  for 
evaluating  the  losses,  are  taken  from  a  WEST  COAST  RESEARCH 
CO.  drawing  [Refs.  7,8]. 

In  each  section  the  loss  K  will  be  a  drop  in  static 
pressure  divided  by  the  dynamic  pressure  at  that  section,  as 
given  by  RAE  &  POPE  [Ref.  2:  p.87]  and  rewritten  below 
Ko  =  K(q/qo), 

K  =  (pi  -  pf  )  /  q  loss  at  each  section, 

Ko  =  coefficient  of  loss  re'^ted  to  test  section, 
pi  =  inlet  section  pressure, 
pf  =  outlet  section  pressure, 
q  =  dynamic  pressure  at  section, 
qo  =  dynamic  pressure  at  test  section. 

In  terms  of  diameter  we  get 

Ko  =  K(Do/D),  where 

Do  =  test  section  diameter, 

D  =  local  tunnel  dicmeter, 

and  finally,  for  energy  ratio  of  the  tunnel, 

ERt  =  jet  energy/summation  of  circuit  losses. 

The  losses  are  evaluated  by  means  of  a  computer  program 
written  in  FORTRAN  IV  language,  which  is  presented  in 
Appendix  A,  and  named  LOSS. 

When  the  evaluation  of  skin  friction  is  needed,  it  is 
made  by  subroutine  FRIC  in  the  LOSS  program,  using  the 
equation  from  RAE  &  POPE  [Ref.  2:  p.88],  shown  immediately 
below. 


1/X 


=  2  logio  ^  ■  ^-8, 


(eqn  2.1) 


where  X 
number . 


skin  friction  coefficient  and  Rey 


Reynolds 


The  Reynolds  number  at  each  section  is  related  to 
Reynolds  number  at  the  test  section,  as  follows: 

Rey  =  VD/j;  . 

From  the  continuity  equation  AV  =  Ao  Vo,  and  hence 
Rey  =  (Do/i;)Vo(Do/D) ,  from  which  is  found 

Rey  =  2.63E  04  Vo(Do/D)  .  (eqn  2.2) 


B.  CYLINDRICAL  SECTIONS 
1 .  Test  Section 

The  losses  for  this  octagonal  section  are  evaluated 
through  the  eqn  2.3,  (eqn  2.44,  Ref.  2,  p.88)  for  the 

equivalent  cylindrical  section, 

Ko=  \(L  /De)(Do/De)^,  (eqn  2.3) 

where  De  is  the  diameter  of  the  circle  whose  circumference 
equals  the  perimeter  of  the  octagon. 

Then,  the  dimensions  for  the  test  section  are: 

Length  =  8.00  feet. 

Diameter  =  4.36  feet. 

This  diameter  is  the  test  section  diameter  (Do)  used 
in  all  the  equations  related  to  the  test  section. 

The  losses  for  the  test  section  at  each  velocity  are 
calculated  by  the  computer  program  LOSS  using  the  equations 
below,  which  were  obtained  from  eqn  2. 3, from  the  equation  of 
the  Reynolds  number,  and  from  the  dimensions  above. 


section  1:  Ko=  1.834X  ,  Rey=  2.63E  04  Vo  . 


2 .  Fan  Duct  Losses 

This  section  begins  after  the  first  corner  with  a 
octagonal  cross  section,  changes  to  a  cylindrical  section 
where  the  two  stages  operate  and  ends  in  another  'octagonal 
section.  From  an  original  wind  tunnel  drawing  (Ref.  8)  we 
have  the  dimensions  for  this  section  and,  using  the  same 
procedure  as  above  (test  section),  we  get  the  following 
dimensions : 

De  =  7.50  ft  L=  9.50  ft  . 

The  same  equations  used  for  test  section  are  used 
here  and  the  equations  for  the  computer  program  are: 

section  2:  Ko=  .  14'48X  ,  Rey=  1.54E  04  Vo  . 

3 .  Constant  Area  Ducts 

The  losses  are  evaluated  in  the  same  manner  as  for 
the  test  section  .  The  dimensions  for  these  two  ducts  are: 

section  3:  De=  13.37  ft,  L=  6.00  ft. 

section  4:  De=  12.73  ft,  L=  6.34  ft. 

The  equations  for  the  computer  program  are: 

section  3:  Ko=  .0051\  ,  Rey=  .86E  04  Vo. 

section  4:  Ko=  .0067X  ,  Rey=  .90E  04  Vo. 

C.  DIFFUSERS 

The  diffusers  are  broken  down  into  parts,  the  first 
(between  test  section  and  first  corner)  into  four,  and  the 
second  (between  the  second  and  third  corners)  into  three, 
the  purpose  of  which  is  to  make  the  evaluation  of  losses 
more  accurately.  The  diffuser  divergence  angle,  which 
strongly  affects  the  evaluation,  is  taken  for  each  part  as 
the  difference  between  the  equivalent  small  and  large 
diameters,  divided  by  the  length. 

The  equation  used  comes  from  RAE  &  POPS  [Ref.  2:  p.89], 
and  it  is 


Ko=  A  (X/8tan(a/2)  + . 6tan(a/2) ) , 


(eqn  2.4) 


where  A  =  ( 1- (D1/D2 )^) (Do/Dl)^  ,  and 

\  =  skin  friction  coefficient  for  Reynolds 
number  given  by  eqn  2.2,  and  based 
on  the  mean  diameter. 


a  =  divergence  angle  between  walls. 
Do  =  test  section  diameter, 

D1  =  smaller  diameter  of  diffuser, 

D2  =  larger  diameter  of  diffuser. 


1 .  Fifst  Diffuser 

The  diffuser  cross  section  from  the  test  section  to 
the  first  corner  varies  along  its  20  feet  of  length,  and  its 
height  H  and  width  W  are  given  as  a  function  of  the  distance 
X  (ft),  along  the  flow  axis,  by  the  equations; 


H  =  3.52  +  .00440  (ft)  , 


(eqn  2.5) 


W  =  5.02  +  .00036  X^  (ft)  .  (eqn  2.6) 

Furthermore,  its  cross-sectional  area  is: 

Area=  .823  W  H  (ft^  )  . 

The  following  table  gives  the  local  values  of  H,  VI, 
Area  and  De  as  function  of  X  for  this  diffuser,  where 

De  =  equivalent  local  section  diameter  (ft). 

The  evaluation  of  the  diffuser  angle  is  made  by  the 
following  equations: 

-For  the  divergence  angle  of  the  equivalent  conical 
diffuser 


-For  the  divergence  angle  between  the  walls  in  the 


TABLE 

I 

PARAMETERS  FOR 

DIFFUSER  1 

X 

H 

W 

Area 

De 

0 

3.52 

5.02 

14.63 

4.32 

5.0 

3.73 

5.06 

15.64 

4.46 

10.0 

4.62 

5.38 

20.60 

5.12 

15.0 

6.44 

6.28 

33.27 

6.51 

18.7 

8.50 

7.38 

51.94 

8.13 

a  =2arc  tan( (D2-Dl)/2L)  . 
vertical  and  horizontal  planes  respectively; 
=arc  tan(  (h2-hl)/2L) , 


(eqn  2.7) 


(eqn  2.8) 


=2arc  tan  ( (w2-wl )  / 2L)  , 


(eqn  2.9) 


where  indexes  1  and  2  refer  to  small  and  large  local 
sections  and  L  is  the  length  between  them. 

For  the  first  diffuser,  the  equations  2.8  and  2.9 
are  used  because  the  angles  between  the  walls  are  different, 
as  shown  in  Table  II  . 

For  a  conservative  analysis,  the  larger  values  are 
chosen  in  computing  the  loss  coefficient  Ko. 

With  the  dimensions  from  this  table  and  the 
equations  2.4  and  2.2  we  got  the  following  equations  to  use 
with  the  computer  program  LOSS  : 


Part 

1 


TABLE  II 

PARAMETERS  FOR  EACH  PART  OF  DIFFUSER  1 


L(ft) 

Dl(ft) 

D2(ft) 

5.0 

4.36 

4.46 

0.52 

2.29 

5.0 

4.46 

5.12 

3.60 

10.23 

5.0 

5.12 

6.51 

9.77 

20.63 

3.7 

6.51 

8.13 

17.50 

30.85 

section 

5: 

Ko=.5367X  + 

.00103, 

Rey=  2.60E 

04 

Vo 

section 

6: 

Ko=.5416A  + 

.0208  , 

Rey=  2.40E 

04 

Vo 

section 

7: 

Ko=.2232\  + 

.0355  , 

Rey=  1.97E 

04 

Vo 

section 

8: 

Ko=.0540A  + 

.0196  , 

Rey=  1.57E 

04 

Vo 

2 .  Second  Diffuser 

For  the  second  diffuser,  the  values  of  D1  and  D2  are 
substituted  into  equation  2.7  for  evaluation  of  the 
divergence  angle,  because  the  angles  between  the  walls  in 
the  vertical  and  horizontal  planes  have  no  large  difference. 


With  the  dimensions  on  Table  III  and  equations  2.4 
and  2.2,  we  got  the  following  equations  to  use  with  the 
computer  program  LOSS. 
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section  9  ;  Ko=.14A  +  .0004 
section  10:  Ko=.09X  +  .0011 
section  11:  Ko=.05\  +  .0030 

D.  CONTRACTION  CONE 


Rey=  1.50E  04  Vo. 
Rey=  1.37E  04  Vo. 
Rey=  l.OOE  04  Vo. 


This  section  takes  the  number  12;  it  is  located  between 
corner  four  and  the  test  section.  The  equation  used  is  from 
RAE  &  POPE  [Ref.  2:  p.91],  and  is 


Ko=  .32\  Lc/Do, 


(eqn  2.10) 


where  \  =  skin  friction  coefficient  for  the  Reynolds  number 
given  by  eqn  2.2  with  D=(Dl+D2)/2;  Lc  =  length  of 
contraction  cone. 

The  equations  for  the  computer  program  are: 
section  12:  Ko=.7338X  ,  Rey=  1.34E  04  Vo. 

E.  SPINNER 

This  section  takes  the  number  13,  and  it  is  located  in 
the  fan  duct.  The  dimensions  are  shown  in  fig  2.2  and  the 
loss  is  evaluated  using  a  method  given  by  Nicolai  [Ref.  10: 
p.  11-21].  The  following  equations  are  for  use  with  the 
computer  program: 

Rel  =  2.45E-04  Vo, 

Cdf  =  .455/(log3^0R«l)’ 

Cdo  =  Cdf  +  .000616  /  Cdf, 

Ko  =  .00837  Cdo, 

Cdo  =  drag  coefficient, 

Cdf  =  friction  drag  coefficient, 

where  Rel  =  Reynolds  Number  related  to  body  diameter. 

The  dimensions  are  : 

L  =  Length  =  12.0  ft, 

D  =  Body  diameter  =  2.0  ft, 

Db  =  Base  body  diameter  =  0.554  ft. 
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Figure  2.2  Spinner  Dimensions. 


F.  CORNERS 

The  corners  are  evaluated  in  three  different  ways: 

(i)  RAE  &  POPE  [Ref.  2:  pp.89,90]. 

(ii)  Bradshaw  and  Pankhurst  [Ref.  5:  p.29  ]. 

(iii)  Pankhurst  and  Holder  [Ref.  3:  p.90]. 

By  RAE  &  POPE  the  loss  equation  is 

Ko=(0.1+(4.55/(log^0  Rn)2-58)(j)o/D)4  ^  2.11) 

where  Rn  =  Reynolds  number  based  on  chord  vane, 

D  =  equivalent  diameter  of  inlet  cross  section. 

The  dimensions  are  in  Table  IV. 

With  the  dimensions  from  Table  IV  and  the  eqn  2.11,  the 
equations  to  be  used  in  the  computer  program  are  found  and 
are  as  follow: 


section 

14: 

Ko  = 

.00690 

+.3122/(logio 

Rn)2 . 58 

Rey  = 

1.646E 

03  Vo  . 

sect  ion 

15  : 

Ko  = 

.00905+ 

.4118/  (log|^Q 

Rn)2-58 

Rey  = 

1.890E 

03  Vo  . 

section 

16; 

Ko  = 

.  00094  + 

. 0428/ ( logjg 

Rn)2-58 

Rey  = 

.  609E 

03  Vo  . 

section 

17  : 

Ko  = 

.00113* 

.0515/(log^0 

Rn)2-38 

Rey  =  .334E  03  Vo. 


section  17: 


jfc.  T  ».  1  WK..  - 


TABLE  IV 

PARAMETERS  OF  CORNERS 

Section 

Chord(ft ) 

D(ft) 

14 

1.042 

8.52 

15 

1.042 

7.95 

16 

1.042 

14.00 

17 

0.521 

13.37 

From  Bradshaw  and  Pankhurst ,  an  empirical  equation  based 
on  tests  at  Reynolds  number  2.0E+05  and  1.9E+06,  gives  the 
loss  coefficient  as 

Cp  =  1.2  (Uc/v)'-^^  ,  (eqn  2.12) 


where  U  =  local  flow  velocity, 
c  =  vane  chord, 

V-  kinematic  viscosity. 

This  equation,  when  related  to  the  test  section  velocity 
by  means  of  equation  Ko=  Cp  (Do/D)^  yields  the  following 
equations  for  use  with  the  computer  program. 


gives  the  variation  of  loss  coefficient  with  the  gap/chord 
ratio  of  the  vanes.  The  tunnel  we  are  dealing  with  has 
corner  -  vane  sections  similar  to  type  b  of  figure  40  of 
reference  3.  Taking  the  dimensions  of  the  corners  and  vanes 
from  the  drawing,  evaluating  the  gap/chord  ratio  for  each 
corner  and  entering  the  graph  on  page  93  of  reference  3  , 
the  value  of  K  (  C  is  used  in  Pankhurst  &  Holder)  is  found. 


section  14: 

Ko  = 

.0130 

Vo" 

,25 

> 

section  15: 

Ko  = 

.0160 

Vo" 

.25 

9 

_ _ 4 

section  16 : 

Ko  = 

.0023 

Vo" 

.25 

9 

section  17 : 

Ko  = 

.0032 

Vo" 

CNJ 

Pankhurst 

and  Holder's 

method  is  a  graphic  one,  which 

■ !  j 
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and  used  with  the  equations  below: 


K=  2  V  H/  p  ,  (eqn  2.13) 

Ko=K  (Do/D)^,  (eqn  2.14) 

where  U  =  local  velocity  at  entry  of  corner, 

VH  =  drop  of  pressure  across  the  corner, 

D  =  equivalent  corner  diameter  at  entry, 

K  =  loss  coefficient  from  the  graph. 

Sample  calculations 

for  corner  i  we  have 
gap  /  chord  =  (8.36  -  2.48)/  12.5  =  0.4704, 

(Do/D)^  =  0.0967, 

K  =  0.235, 

Ko  =  0.0227. 

The  losses  in  all  four  corners,  as  given  by  Pankhurst 
and  Holder's  method  are  shown  in  Table  V  . 

G.  SUMMARY  OF  LOSSES 

The  losses  evaluated  by  the  computer  program  for 

sections  1  through  17  for  a  test  section  velocity  of  200 
Knots,  are  shown  in  Table  VI  ;  values  for  velocities  from 

100  to  200  Knots  are  shown  in  Table  VII  . 

The  energy  ratio  given  by  Ert  is  multiplied  by  .9  in 
order  to  take  into  account  the  losses  due  to  leaks  and 
joints . 

The  results  presented  in  Table  VII  include,  for  purposes 
of  comparison,  the  losses  evaluated  by  the  other  two 

methods.  At  the  design  airspeed  of  200  Knots,  Pankhurst  & 
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TABLE  V 

CORNER  LOSSES  BY  PANKHURST  AND  HOLDER 


corner 

gap/ chord 

(Do/D)^ 

K 

Ko 

1 

0.4704 

0.0967 

0.23 

0.023 

2 

0.4256 

0.1334 

0.27 

0.036 

3 

0.4752 

0.0094 

0.23 

0.002 

0.4656 

0.0113 

0.24 

0.003 

Holder's  method  predicts  that  245  HP  will  be  required,  while 
Bradshaw  &  Pankhurst  (a  later  work)  indicates  a  need  for  170 
HP.  Rae  &  Pope's  method  requires  196  HP.  The  power  input  to 
the  tunnel  is  300  x  fan  eff iciency( . 85^)  =  217  HP. 

The  description  of  the  computer  program  LOSS  is 
presented  in  Appendix  D. 


TABLE  VI 

LOSS  FOR  EACH  SECTION  OF  TUNNEL  AT  V0=200  KNOTS 


SEC 

ReylO"® 

X 

Ko 

%Loss 

1  test 

8.88 

0.0082 

0.0151 

10.99 

2  cyl 

5.16 

0.0090 

0.0013 

0.95 

3  cyl 

2.90 

0.0100 

0.0001 

0.04 

4  cyl 

3.04 

0.0100 

0.0001 

0.05 

5  dif 

8.77 

0.0083 

0.0055 

3.99 

6  dif 

8.07 

0.0084 

0.0253 

18.42 

7  dif 

6.65 

0.0090 

0.0375 

27.25 

8  dif 

5.29 

0.0090 

0.0201 

14.59 

9  dif 

5.06 

0.0090 

0.0017 

1.21 

10  dif 

4.62 

0.0091 

0.0019 

1.39 

11  dif 

3 . 37 

0.0100 

0.0035 

2.58 

12  con 

4.53 

0.0091 

0.0067 

4.86 

13  spi 

0.0008 

0.60 

14  cor 

4.56 

0.0067 

4.88 

15  cor 

4.86 

0.0092 

6.70 

16  cor 

2.77 

0.0009 

0.69 

17  cor 

2.90 

0.0011 

0.82 

120.0  38.0835  43.5826  53.7830 
140.0  60.1842  68.9301  85.2725 
160.0  89.2441  102.3171  126.8912 

143.8930  179.1000 
195.9971  244.5512 


180.0  125.2569 

200.0  170.4057 


III.  BLADE  PERFORMANCE 


A.  INTRODUCTION 

The  "academic"  wind  tunnel  has  two  stage  counter  - 
rotating  fans,  each  fed  by  an  electric  motor  of  150  hp  and 
with  a  fixed  speed  of  1200  rpm.  The  main  steps  used  in  the 
design  procedure  were  found  in  the  Larson  report  [Ref.  4,] 
which  contains  the  original  blade  design.  The  equations  and 
the  assumptions  made  during  the  design  indicate  that  three 
dimensional  effects  were  not  accounted  for.  For  this  reason 
the  method  given  by  Borst,  which  is  based  on  the  blade 
element  approach  and  cascade  theory  for  determining  these 
three  dimensional  effects,  will  be  used  to  predict  the 
performance  of  the  blades.  Usefully,  this  method  requires 
data  of  two  dimensional  airfoils  only,  in  order  to  determine 
the  force  on  each  section  of  the  blade. 

So,  the  blade  system  with  all  assumptions  made  during 
the  design,  (see  section  B  -  Larson  report),  and  the  Borst 
method,  (see  section  C  ),  are  the  material  necessary  to 
evaluate  the  performance  of  the  fans.  Knowing  the  lift 
developed  at  each  blade  station  (eqn.  3.11)  one  can  find  the 
torque  (eqn.  3.12)  at  each  station.  A  blade  station  is 
defined  non-dimensionally  by  the  ratio  of  the  radius  of  that 
station  to  the  blade  tip  radius;  the  root  station  for  this 
tunnel  is  at  .267  and  the  tip  station  is  at  unity.  Then 
adding  the  torques  at  all  stations, and  multiplying  by  the 
rotational  speed,  one  gets  the  power  for  the  blade. 

B.  THE  LARSON  REPORT 

This  report  describes  the  procedure  that  was  followed  to 
design  the  blades,  together  with  the  assumptions  made  during 
the  calculations. 

The  most  important  assumptions  were: 

-  85%  fan  efficiency  for  each  stage. 


Axial  velocity  distribution  of  the  flow 
entering  the  first  stage  fan  is 

Cm= . 12  Cmr  (rt-r)+Cmt  ,  (eqn  3.1) 

where  Cmr  =  axial  velocity  at  blade  root, 

Cmt  =  axial  velocity  at  blade  tip, 
rt  =  tip  radius, 
r  =  local  radius. 

-  Maximum  test  section  velocity  equal  to  200 
knots . 

-  Design  lift  coefficient  equal  to  0.64,  leading 
to  the  requirement  that  Cl=  .8  . 

-  Airfoil  type  NACA  16  -  X12  . 

This  design  was  based  on  forced  vortex  flow,  in  other 
words,  the  difference  in  static  pressure  across  the  vortex 
is  equal  to  the  difference  in  tangential  velocity  (dynamic) 
pressure.  This  leads  us  to  the  loading  parameter  (solidity  x 
lift  coefficient)  and  the  mean  velocity  angle,  as  one  can 
see  in  Appendix  A,  whose  equations  are: 

a  Cl  =(  2(A-B)C  )/(  l+(Cd/Cl)  D,  and  (eqn  3.2) 


Pm  =(;3i+/32)/2  ,  where  (eqn  3.3) 

A=  tan)3]^  , 

B=  tanj32  ’ 

C=  cosj3n,  , 

D-  tan^j^  . 

Using  these  equations  and  the  assumptions  mentioned 
above,  the  dimensions  of  the  blades  were  calculated,  and  are 
given  on  the  drawing  of  the  blades  (Ref.  9).  The  blades  were 
divided  into  11  (eleven)  stations  from  root  to  tip;  for  each 
station  the  chord,  thickness,  leading  edge  angle,  etc... 


were  determined. 


C.  THE  BORST  METHOD 

This  method  is  based  on  the  blade  element  approach  and 
the  vortex  theory  for  determining  the  three-dimensional 
effects,  so  that  two-dimensional  airfoil  data  can  be  used 
for  determining  the  resultant  force  on  each  blade  element. 
By  the  momentum  flow  theory,  Borst  obtained  the  loading 
parameter  (solidity  x  lift  coefficient),  as  a  function  of 
the  angles  j3]^,  ^2  /^m  ’  that  are  dependent  upon  an 

induced  angle  of  attack  (a^)  at  values  of  constant  inflow 
angle . 

The  loading  parameter  is: 


ffCl=(  2  B  (C  -D)  )  /  E, 


.*  v^- 


5&V5: 

.V.--  .N 


(eqn  3.4) 


W-*  - 

..  .J 

V  ^  • 


where  A  =  cos^^  , 

B  =  cosy  , 

C  =  tanjS^  , 

D  =  tanj32  > 

E  =  cos(j3^-y), 

Cl  =  lift  coefficient, 

Cd  =  drag  coefficient, 

=  apparent  inlet  angle, 

^2  -  apparent  outlet  angle, 

=  angle  of  the  mean  velocity  vector, 
y  =  atan  (Cd/Cl), 
a  -  solidity. 

Now,  based  on  the  equations 


^  • 
f  •  • 


( eqn  3.5) 


^m  “  ^  1  '  “  i 


( eqn  3.6) 


V 


and  the  two  dimensional  airfoil  data  and  design  parameters 
for  the  blade,  one  can  iterate  the  equation  3.4  and  get  the 
induced  angle  of  attack  for  a  specified  test  section 
velocity  and  blade  station. 

The  angle  of  attack  of  the  blade  at  any  given  station 


a  =  -j3  +  fl  , 


(eqn  3.7) 


where  /3=  stagger  angle  (Ref.  9,  Ref.  4  ,  p.7),  which  is  the 
angle  between  the  chord  line  of  the  blade  station  and  the 
line  parallel  to  the  axis  of  rotation  of  blades,  (see  fig. 
3.1  )  and 

8=  pitch  angle  (adjustable). 

The  Figure  3.1  shows  the  relation  between  the  variables 
given  above. 

D.  TORQUE  AND  POWER  EVALUATION 

We  will  use  the  Borst  '  method  (Ref.  1),  to  get  the 
induced  angle  for  correction  of  the  two  dimensional  angle  of 
attack,  and  the  Larson  Report  (Ref. 4)  to  obtain  all  the 
dimensions  of  the  blades  and  the  preliminary  design 
conditions . 

With  this  material,  we  are  able  to  follow  the  procedure 
outlined  below  to  get  to  the  value  of  torque  and  power  for 
each  fan  stage.  This  procedure  will  lead  us  to  all 
velocities  and  angles  of  each  blade,  finalizing  with  the 
value  of  the  torque  on  each  blade  section.  After  that,  we 
evaluate  the  torque  for  the  whole  blade,  and  finally  the 
torque  and  power  for  the  fan  stage. 

PROCEDURE  FOR  EVALUATION  OF  TORQUE  AND  POWER 

STEP  1)  From  eqn.  3.1  for  a  given  test  section  velocity 
and  n=  1,200  rpm,  calculate  W1  =(Cm^  +  U^  )l/2, 
where  U=  2^  r  n. 

As  one  can  see,  W1  =  f(  r,Vo)  since  Cm  =  f(  i,Vo).  In 
Appendix  A  the  axial  flow  velocity  (Cm)  was  developed  as 


Figure  3.1  Angle  of  .Attack  Relations. 


function  of  Vo. 

STEP  2)  Calculate  =  arc  tan  (  U  /  Cm). 

STEP  3)  From  the  drawing  of  the  blades  (Ref.  9),  the 
values  of  chord  are  found  and  then  the  solidity  for  the 
blade  section  is  obtained  from 

a  =  B  c  /  2  TT  r  ,  (eqn  3.8) 

where  B  =  number  of  blades,  c  =  station  chord  and  r  = 
station  radius . 

STEP  4)  From  the  "Aerodynamic  Characteristics  of  24  N’.ACA 
16  -  Series  Airfoils"  -  TN  1546  ,  (Ref. 12,  p.  49,  50)  the 
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lift  and  drag  coefficient  for  each  section  are  available  at 
MACH  number  M  =  W1  /  a  ,  where  a  =  sound  velocity  at 
temperature  of  100  degrees  Fahrenheit,  and  for  the  design 
lift  coefficient  of  0.64,  which  was  assumed  by  Larson  (Ref. 
4  ,p,5,6).  The  angle  of  attack  used  to  enter  the  Cl  vsa 
curve  is  found  from  eqn.3.7.  Find  Cl,  Cd  and  evaluate  y  = 
arc  tan(Cd/Cl). 

STEP  5)  Assuming  an  initial  value  of  induced  angle, 
evaluate  the  angles  given  by  equations  3.5  and  3.6. 

STEP  6)  Evaluate  the  right  side  of  equation  3.4  using 
the  values  found  in  steps  1,2,  and  5;  evaluate  the  left  side 
of  equation  3.4  using  the  values  found  in  steps  3  and  4.  If 
both  sides  are  equal,  the  value  of  induced  angle  assumed  is 
correct,  otherwise,  iterate. 

STEP  7)  The  corrected  angle  of  attack,  found  by  using 
the  induced  angle,  is  seen  in  figure  3.1  and  evaluated  as 
follows ; 

a^t.  =  a  -  .  (eqn  3.9) 

STEP  8)  With  corrected  angle  of  attack,  recalculate  the 
step  4  above  and  get  the  corrected  lift  and  drag 
coefficient,  and  evaluate  y  with  these  values. 

STEP  9)  The  mean  velocity  vector  and  the  lift  for  each 
blade  station  are  evaluated  using  the  following  equations; 

Wm  =  Cm  /cos/3^  ,  (eqn  3.10) 


Li  =  1/2  p  Wm^  Cl  c  , 


(eqn  3.11) 


where  Li  =  lift  at  blade  station  i  per  unit  length. 
Cl  =  lift  coefficient  at  blade  station, 

Wm  =  mean  velocity  vector  at  blade  station, 
c  =  chord  at  blade  station. 


.-v 
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STEP  10)  By  definition,  the  lift  force  is  perpendicular 
to  Wm.  For  evaluation  of  the  torque  it  is  necessary  to  get 
the  component  of  the  resultant  force  in  the  plane  of  the 
fan,  as  shown  in  figure  3.2.  Then, 

Ri=  Li/ cos.  y  , 

Fi  =  Ri  cos  (j3nj  -y  ), 

where  y  =  arc  tan  (Cd/Cl).  With  Cd  and  Cl  corrected  and 
evaluated  at  step  8),  and  Li,Ri,Fi  are  respectively  the  lift 
force,  the  resultant  force  of  the  lift  and  drag,  and  the 
force  in  the  plane  of  the  fan. 


Figure  3.2  Lift  Developed  at  Each  Blade  Station. 

STEP  ll)The  force  per  unit  length,  Fi,  averaged  between 
the  value  at  one  radial  station  and  next,  times  the 
difference  in  radius  between  these  stations,  is  the  force 
used  for  torque  evaluation.  This  force  times  the  means 
radius  of  these  two  stations  is  the  torque.  At  this  point  it 
is  necessary  to  go  through  step  1  to  10  again,  in  order  to 
get  Fi  for  the  next  station.  Then  the  torque  on  a  blade 
section  between  two  consecutive  stations  is: 


Ti  +  Fi  )  S  R  /  2 


(eqn  3.12) 


where  S  =  -  r^)  length  between  stations, 

R  =  +  T^)/2, 

and  i  refers  to  the  blade  station. 

STEP  12)  Adding  the  section  torques  along  the  blade, 
multiplying  the  result  by  the  number  of  blades  and  by  the 
rotational  speed,  the  power  required  is  found.  This  power 
provides  the  test  section  velocity  which  was  specified  at 
the  beginning  of  the  procedure. 

E.  INTERACTION  BETWEEN  STAGES 

For  evaluating  the  torque  and  power  for  the  second 
stage,  the  procedure  given  in  the  previous  section  can  be 
used,  but  with  some  adjustments  for  the  different  geometric 
characteristics  of  the  blade;  for  example,  the  chord, 
stagger  angle,  etc.  .  The  tangential  velocity  used  is  the 
tangential  velocity  of  the  blade  at  the  given  station  plus 
the  tangential  velocity  of  the  outflow  of  the  first  stage  at 
that  station.  It  should  be  noted  that  the  axial  inflow 
velocity  is  assumed  the  same  as  for  the  first  stage.  In 
figure  3.4  this  interaction  can  be  seen. 

The  second  stage  should  be  able  to  accept  the  tangential 
velocity  introduced  in  the  flow  by  the  first  stage  and  add 
the  same  amount  of  power  as  added  by  the  first  stage.  The 
flow  leaving  the  second  stage  then  has  no  rotation. 


for  Torque  Evaluation  at  Each  Station. 


Figure  3.4  Axial  Flow  Veloc 
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F.  CONCLUSION 


Using  the  computer  program  BORST  and  the  procedure 
discussed  in  this  chapter,  Tables  VIII  through  XX  were 
obtained.  For  all  these  tables,  a  pitch  angle  of  2.4 
degrees  was  used. 

The  Table  VIII  shows  the  values  of  the  velocities  from 
root  to  tip,  for  Vo=  200  Knots,  of  the  first  fan  stage. 
These  velocities  are  defined  as  follows: 

U  =  tangential  blade  velocity  (ft/sec). 

Cm  =  axial  flow  velocity  (ft/sec), 

W1  =  apparent  flow  velocity  (ft/sec). 

The  chord  at  each  station  is  given  in  feet. 

The  Table  IX  shows  the  Reynolds  and  Mach  numbers  at  each 
station,  for  the  first  fan  stage  and  at  V0=  200  Knots.  It 
is  clear  that  compressibility  effects  appear  at  the  outer 
20%  of  the  blade  only.  The  maximum  Mach  number  was  less  than 
.45  and  the  values  of  Cl  vs  Alpha  and  Cl  vs  Cd  were  obt^ained 
by  linear  interpolation  between  the  graphs  for  M= . 3  and 
M=  .45. 

The  Table  X  shows  the  angles  between  the  axial  flow 
velocity  vector  and  the  velocities  relative  to  the  blade, 
the  stagger  angle,  the  induced  angle,  the  angle  of  attack, 
and  the  lift  and  drag  coefficients.  These  values  are  for  the 
first  fan  stage,  at  each  station  and  at  Vo=  200  Knots.  The 
values  of  lift  coefficient,  are  distributed  along  the 
stations  of  the  blade  from  .662  to  .896,  around  the  uniform 
value  of  .8,  assumed  by  Larson.  These  variables  are  defined 
in  Figures  3.1,  3.4,  and  the  angles  are  given  in  degrees.  As 
expected,  at  the  root,  the  induced  angle  is  the  largest, 
because  the  greatest  interference  between  the  blades  occurs 
there . 

The  Table  XI  shows  the  angle  of  attack,  lift  and  drag 
coefficients,  all  with  correction  for  the  induced  angle. 
Also  shown  are  the  lift,  the  resultant  force  and  the  torque. 
This  Table  refers  to  each  station  of  the  first  fan  stage, 
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with  Vo=  200  Knots.  The  values  of  corrected  Cl  (CLCR)  are 
distributed  along  the  blade  from  .358  to  .709,  differently 
from  that  assumed  by  Larson,  which  was  a  uniform  .8  . 

The  variables  in  Table  XI  are  defined  as  follow: 

ALPCR  =  corrected  angle  of  attack  (degrees), 

CLCR  =  corrected  lift  coefficient, 

CDCR  =  corrected  drag  coefficient, 

DELL  =  lift  force  per  unit  length  (Ibf/ft), 

DELF  =  tangential  force  per  unit  length  (Ibf/ft), 

DTOR  =  torque  at  middle  of  two  stations  (Ibf-ft). 

The  Tables  XII,  XIII  and  XIV,  apply  to  the  second  fan 
stage.  The  variables  have  the  same  definitions  as  for  the 
first  stage.  The  Table  XII  shows  the  inflow  velocity 
relative  to  the  blade  (W3)  larger  than  Wl.  This  is  due  to 
the  tangential  flow  velocity,  which  has  a  finite  value  after 
the  first  fan  stage.  At  the  root  the  difference  between  W3 
and  Wl  is  larger  than  at  the  tip,  this  means  that  the 
tangential  flow  velocity  at  the  root  is  larger  than  at  the 
tip.  The  Table  XIII  shows  the  Mach  number  less  than  .45,  and 
the  interpolation  described  above  was  again  used.  The  Table 
XIV  shov?s  the  Cl  varying  around  .8,  that  is,  from  .875  to 
.710  as  for  the  first  fan  stage. 

The  Table  XV  refers  to  the  second  stage,  with  the  same 
variables  as  Table  XI  plus  the  variable  BETA1-BETA4.  This 
Table  shows  the  difference  between  BETAl  and  BETA4 ,  which  is 
around  zero.  This  means  that  the  flow  after  the  second  fan 
stage  has  no  tangential  velocity. 

The  Tables  XVI  and  XVII,  are  for  the  uncorrected  values 
of  the  angle  of  attack,  for  first  stage  at  Vo=200  Knots.  The 
Table  XVI  shows  the  values  of  Cl  around  .8  .  The  Table  XVII 
shows  values  of  torque  (DTOR)  larger  than  those  with  the 
corrected  angle  of  attack  (Table  XI). 

The  Tables  XVIII,  XIX  and  XX  are  for  the  second  fan 
stage  without  correction  of  the  angle  of  attack.  A 
comparison  of  the  velocities  W3  from  Table  XVIII,  and  W3 


from  Table  XII  shows  that  the  difference,  due  to  the 
correction  of  the  angle  of  attack  has  affected  the  outflow 
of  first  fan  stage.  The  Table  XIX  shows  the  variation  of 
BETAS  as  a  result  of  the  variation  of  W3 ,  just  mentioned 
above.  The  variation  of  Cl  was  not  so  large  as  to  affect 
the  value  of  Cd.  A  comparison  of  Tables  XV  and  XX  shows 
that  the  interference  between  the  blades  reduces  the  torque 
by  about  1/3  at  the  root  and  very  little  at  the  tip.  This  is 
further  discussed  below. 

In  order  to  visualize  the  effects  of  the  correction  on 
the  angle  of  attack,  two  graphs,  each  with  four  curves,  were 
plotted;  these  are; 

1  Figure  3.5;  the  tangential  force  on  the  blade  is  shown 
for  each  station,  at  a  specified  test  section  velocity 
of  200  Knots  and  at  a  pitch  angle  of  2.4  degrees.  The 
curves  are  referred  to  the  first  and  second  stages  of 
the  blades,  with  and  without  correction  of  the  angle  of 
attack.  This  figure  shows  again  the  effects  of  the 
interference.  The  tangential  forces  have  the  largest 
difference  at  the  root,  and  are  practically  equal  at 
the  tip.  The  data  were  taken  from  Tables  XI  ,  XV  , 
XVII  and  XX  . 

2  Figure  3.6;  the  angle  of  attack  is  shown  for  each 
station  of  the  blade,  at  a  specified  test  section 
velocity  of  200  Knots  and  at  a  pitch  angle  of  2.4 
degrees.  The  corrected  angle  of  attack  at  the  root 
became  negative  and,  at  the  tip,  that  angle  is 
practically  equal  to  the  uncorrected  one.  This  again 
shows  the  strong  effects  of  interference  at  the  blade 
root.  The  curves  are  referred  as  above,  and  the  data 
were  taken  from  Tables  X  ,  XI  ,  XIV  and  XV  . 


Figure  3.5  Tangential  Force  on  the  Blade  at  Each  Station 
Vo=200  Knots,  Pitch  Angle=  2.4  degrees. 


TABLE  VIII 


VELOCITIES  AT  EACH  STATION, 
FIRST  STAGE,  V0=200  KNOTS 


STATION 

CHORD 

U 

CM 

W1 

0.267 

1.112 

125.821 

152.606 

197.787 

0.300 

1.049 

141.372 

150.340 

206.369 

0.350 

0.958 

164.934 

146.907 

220.873 

0.400 

0.875 

188.496 

143.475 

236.887 

0.450 

0.802 

212.058 

140.042 

254.126 

0.500 

0.738 

235.620 

136.609 

272.357 

0.600 

0.631 

282.744 

129 . 743 

311.090 

0.700 

0.546 

329.868 

122.877 

352.010 

0.800 

0.481 

376.992 

116.011 

394.438 

0.900 

0.422 

424.116 

109.145 

437.935 

1.000 

0.344 

471.240 

102.279 

482.212 

TABLE  IX 

REYNOLDS  AND  MACH  NUMBERS  AT  EACH  BLADE  STATION, 
FIRST  STAGE,  V0=200  KNOTS 


STATION 

REY 

MACH 

0.267 

1324694.000 

0.171 

0.300 

1303483.000 

0.178 

0.350 

1274545.000 

0.190 

0.400 

1249081.000 

0.204 

0.450 

1227461.000 

0.219 

0.500 

1210513.000 

0.235 

0.600 

1181581.000 

0.268 

0.700 

1158453.000 

0.303 

0.800 

1143157.000 

0.340 

0.900 

1112776.000 

0.378 

1.000 

998700.900 

0.416 
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TABLE 

XII 

Vr.LOCITTES  AT 
SECOND  STAGE, 

EACH  STATION, 

V0=200  KNOTS 

STATION 

CHORD 

U 

CM 

W3 

0.267 

0.931 

125.821 

152.606 

238.513 

0.300 

0.893 

141.372 

150.340 

242.660 

0.350 

0.834 

164.934 

146.907 

252.430 

0.400 

0.777 

188.496 

143.475 

263.640 

0.450 

0.724 

212.058 

140.042 

276.489 

0.500 

0.673 

235.620 

136.609 

292.340 

0.600 

0.590 

282.744 

129.743 

327.947 

0.700 

0.519 

329.868 

122.877 

366.731 

0.800 

0.461 

376.992 

116.011 

407.820 

0.900 

0.411 

424.116 

109.145 

450.087 

1.000 

0.373 

471.240 

102.279 

492.649 

TABLE  XIII 


REYNOLDS  AND  MACH  NUMBER  AT  EACH  STATION, 
SECOND  STAGE,  VO=  200  KNOTS 


STATION 

REY 

MACH 

0.267 

1337683.000 

0.206 

0.300 

1305394.000 

0.209 

0.350 

1268230.000 

0.218 

0.400 

1234027.000 

0.227 

0.450 

1205890.000 

0.238 

0.500 

1185210.000 

0.252 

0.600 
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ANGLES,  LIFT  AND  DRAG  COEF.  AT  EACH  STATION 
WITHOUT  CORRECTION  OF  ANGLE  OF  ATTACK,  FIRST  STAGE,  V0=200  KNOTS 
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IV.  TUNNEL  PERFORMANCE 


A.  INTRODUCTION 

The  tunnel  performance  depends  on  losses  through  the 
tunnel  circuit  and  on  the  fan  efficiency. 

The  performance  will  be  evaluated  to  ascertain  how  the 
loss  around  the  tunnel  circuit  is  related  to  the  velocity  at 
the  test  section,  and  if  the  fans  can  provide  adequate 
energy  to  generate  this  test  section  velocity.  The  energy 
provided  by  the  blade  system  depends  upon  the  pitch  angle 
which  must  be  limited  in  order  to  avoid  stall.  The  pitch 
angle  is  found  for  the  first  and  the  second  stages  and 
checked  if  the  flow  leaving  the  second  stage  has  zero 
rotational  velocity  as  required. 

The  tunnel  performance  will  be  shown  through  the 
operational  envelope  of  the  -blade  system. 

B.  TUNNEL  LOSSES 

As  we  have  seen  in  CHAPTER  TWO  the  losses  of  each 
section  depend  on  the  flow  velocity,  which  implies  that  they 
depend  upon  of  test  section  velocity 

Relating  the  losses  for  each  tunnel  section  to  the  test 
section  velocity  we  can  get  the  total  losses  as  a  function 
of  this  velocity,  and  then  relate  the  latter  to  flow  energy 
required. 

Using  the  program  LOSS  with  test  section  velocity 
varying  from  100  to  200  Knots,  by  increments  of  20  Knots,  we 
get  the  figure  4.1  which  shows  the  energy  required  (rjbhp .  ) 
versus  test  section  velocity  (Vo)  . 

C.  OPERATIONAL  ENVELOPE 

Following  the  method  given  in  Chapter  3,  we  find  the 
power  required  from  the  fans  as  a  function  of  test  section 
velocity,  with  a  number  of  pitch  angles  of  the  first  and 
second  stage  blades. 


Figure  4.1  Energy  Required  to  Overcome  Tunnel  Losses. 

Since  we  know  the  losses  around  the  tunnel  in  relation 
to  test  section  velocity,  we  can  adjust  the  pitch  angle  of 
blades  until  we  get  the  energy  coming  into  the  flow  to  match 
the  losses  at  same  test  section  velocity.  The  program  BORST 
was  used  for  this. 


In  order  to  test  the  sensitivity  of  the  blade  flow  to 
the  approaching  velocity  profile,  we  compare  the  performance 
of  the  blades  in  uniform  and  skewed  approach  flow. 

Using  the  axial  flow  velocity  profile,  originally  used 
by  Larson  and  given  by  equation  4.1  with  the  program  BURST, 
we  got  the  data  shown  in  Table  XXI  .  This  Table  presents 
the  pitch  angle  for  a  specified  energy  and  test  section 
velocity,  the  maximum  and  minimum  pitch  angle  at  that  test 
section  velocity  and  its  related  energy.  These  values  were 
plotted  and  are  shown  in  Figure  4.2  .  This  Figure  has  three 
limit  curves,  root  stall,  tip  stall  and  fan  efficiency.  The 
root  stall  curve  limit  is  the  locus  of  points  at  which,  at  a 
specified  test  section  velocity  and  pitch  angle,  the  blade 
stalls  at  the  root.  The  tip  stall  curve  limit  is  for  stall 
at  the  tip.  The  fan  efficiency  curve  limit  is  the  locus 
where  the  operation  of  the  fans  is  limited  by  the  input 
power  and  the  (assumed)  fan  efficiency,  that  is,  300  hp  x 
.85^  =  217  hp.  This  means  that  the  power  provided  by  the 
fan  system  can  not  exceed  217  hp  approximately.  The  pitch 
angles  for  both  stages  are  equal  at  the  design  speed  of 
Vo=200  Knots. 

Using  the  uniform  flow  velocity  given  by  equation  4.2  , 
we  found  the  data  shown  in  Table  XXII  .  This  Table  shows 
that  the  fans  can  operate  without  stall,  only  at  a  section 
test  velocity  from  160  to  180  Knots.  These  values  were 
plotted  in  Figure  4.3  .  This  Figure  has  only  the  curves  of 
root  and  tip  stall.  The  fans  stall  before  reaching  the 
maximum  power  allowed  by  the  fan  efficiency.  The  blades  of 
second  fan  stage  are  practically  stalled  at  the  root. 

The  axial  velocity  distribution  specified  for  Table 
XXI  and  Table  XXII  are  developed  in  Appendix  A  and  shown 


Cm  =  .2034  (  1  -  X  )  +  .303  Vo  ,  and 


(eqn  4.1) 


Cm  =  .363  Vo.  (eqn  4.2) 

D.  FINAL  CONCLUSION 

The  large  subsonic  wind  tunnel  at  NFS,  powered  by  two 
counter-rotating  fans,  never  achieved  the  design 
specifications.  The  reason  for  the  poor  performance  of  the 
tunnel  was  unknown  but  believed  to  be  due  to  either  poorly 
designed  fan  blades  or  to  separation  in  the  diffuser.  Being 
easier  to  analyse,  the  fan  blades  were  chosen  for  initial 
study.  However,  to  do  this  analysis,  it  was  necessary  to 
compute  the  losses  in  the  whole  tunnel  circuit. 

The  results  can  be  summarized  as  follows: 

•1)  Most  of  the  losses,  about  60%  of  the  total,  was  in  the 
first  diffuser.  The  exit  end  of  this  diffuser  has  a  high 
value  of  divergence  angle  in  the  vertical  plane,  which  can 
be  the  cause  of  a  possible  flow  separation. 

2)  V/hen  the  inflow  velocity  profile  to  the  fans  is  assumed 
uniform,  the  no-stall  operational  envelope  is  substantially 
reduced,  so  there  is  a  much  greater  likelihood  of  blade 
stall . 

3)  The  variation  of  the  pitch  angle  affects  the 
performance  of  the  fans.  The  pitch  angles  for  both  stages 
should  be  practically  the  same  in  order  to  have  no  flow 
rotation  after  the  outlet  of  the  second  stage,  provided  the 
inflow  velocity  profile  is  the  same  as  that  assumed  in  the 
design. 

4)  In  conclusion,  although  this  new  blade  element  method 
used  to  analyse  the  fan  flow  predicts  some  deviation  from 
the  original  isolated  blade  analysis,  it  was  found  that  the 
original  design  was  adequate  and,  therefore,  the  tunnel 
problem  is  most  probably  due  to  separation  in  the  diffuser. 


:.-Vvv^.r,v-v.v>. 


59 


TABLE  XXI 

OPERATIONAL  ENVELOPE  FOR  EQUATION 
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APPENDIX  A 

SPECIFIED  FUNCTIONAL  FORMS 


CORNER 

The  evaluation  of  the  corner  loss  using  the  turning  vane 
loss  function  presented  in  appendix  B  of  reference  6  , 

(p.35)  with  the  flow  turning  angle  equal  to  90  degrees, 
leads  to  the  equation  given  by  RAE  &  POPE  [Ref.  2:  p.90  ]. 
The  equation  for  the  turning  vane  loss  coefficient  is 

KTV=-1.605670E-01^-1.446753E-02^  -2 . 570748E-04<^2  +  2 . 066207 

E-06<^,^  -6.335764E-O90^  , 

with^=  90  degrees  is  found  KTV=  .15. 

The  equation  for  turning  vane  loss  function  is 
K  =  KTV  (2/3+(l/3(logi^Q  Rnr  /log^Q 
where 

Rnr  =  500,000  (reference  Reynolds  number, 
according  to  ref.  2  p.89  ), 
and  loS]_o  500,000  =  5.7  . 

Substituting  these  values  in  the  equation  of  the  turning 
vane  loss  function,  the  equation  given  by  Rae  &  Pope  appears 
as  follows: 

K  =  0.15  (2/3  +  29.7/logj^0Rn)2-58)^ 

K  =  (0.1  ^  4.45/  ), 

Ko  =  K  (Do/D)^  . 

LOADING  PARAMETER  (Larson) 

The  loading  parameter  used  by  Larson  in  his  report  comes 
from  blade  element  theory,  where  the  effects  o"  drag  are 
considered.  This  is  also  shown  in  Osborne  [Ref.  11:  p.l47], 
by  equation 

Vu/Wm=  1/2  Cl(  1+  Cd/Cl  cotan  (^-a)) 
where  =  mean  blade  velocity  angle  . 

The  Larson  report  (pp.  5,6)  gives 

cotan  (j3  -a)  =  tan  (jSnj) , 
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and  rearranging  the  term  Vu/Wm  as 

Vu/Wm  =  (U/Cm  -  (U-Vm)/Cm)  Cm/Wm  ,  and  using  the  axial 
velocity  triangles,  the  following  relations  are  assured: 
U/Cm=  tanjS^,  (U-Vu)/Cm  =  tanj32  ’ 

Cm/Wm  =  cosjS^i 

Substituting  these  relations  into  the  equation  given  by 
Osborne,  the  final  equation  is 

aCl  =  (2  (tan^j^- tanjS2  )cosjSni)/ ( 1+Cd/Cl 
equation  used  by  Larson. 

AXIAL  FLOW  VELOCITIES. 

Larson  assumed  in  his  report  the  axial  flow  velocity  as: 

Cm  =  .12  rt  (1-  r/rt)  +  Cmt  (1) 

where  rt=  tip  radius,  r=  station  radius,  Cmt=  axial 
velocity  at  tip  blade,  and  Cmr  =  axial  velocity  at  the  root 
of  the  blade.  Integrating  the  product  of  the  cross  section 
area  of  the  fan  duct  at  the  blades  times  Cm,  from  root  to 
tip,  and  setting  it  equal  to  the  product  of  the  area  times 
the  velocity  of  the  test  section,  we  find: 

14.9  Vo  =  5.465  Cmr  +  41.037  Cmt. 

In  the  equation  of  Cm,  when  r  is  the  radius  of  the  blade 
root.  Cm  =  Cmr,  and 

Cmr  =  Cmt  /  0.67  .  (2) 

Equating  (1)  and  (2),  the  axial  velocity  profile  for  our 
tunnel  is  found  as: 

Cm  =  (  0.2034  (  1-  r/rt  )  +  .303  )  Vo 
UNIFORM  AXIAL  FLOW  VELOCITY. 

This  axial  profile  velocity  is  used  only  for  purposes  of 
comparison.  From  the  continuity  equation  Cm  =  Ao  Vo  /  A 
and,  taking  the  known  values,  we  find  Cm  =  .363  Vo 


APPENDIX  B 
SAMPLE  CALCULATIONS 


Evaluation  of  the  blade  station  torque. 

For  first  stage 

Following  the  'torque  and  power'  procedure  given  in 
chapter  three,  at  a  specified  station  the  data  are  as 
follows: 

General  data: 

-  curves  of  Cl  vsa  and  Cl  vs  Cd  for  air¬ 
foil  NACA  16-X12  at  M= . 3  and  M=.45. 

a=l,160  ft/sec 

-  p  -  .0023  slug/ft^ 

-  V  =.000166  ft^/sec 

-  rt  =  3.75  ft 

-  Cldes  -••64 

-  Vo  =  200  Knots  =  337.55  ft/sec 

-  B  =  4  blades 

-  9  =  2.4  deg. 

At  specified  station: 

-  X  =  .3 

-  c  =  1.049  ft 


-  =  40.37  deg. 

STEP  1)  Cm  =  (. 2034(1-. 3)+.303)x337. 55  =  150.34  ft/sec 
U  =  2x3. 1416x. 3x3. 75x1200/60  =141.37  ft/sec 
W1  =  (150.342+141.372)1/2  =  206.37  ft/sec 
STEP  2)  j3i  =  arc  tan(141. 37/150 . 34)  =  43.24  deg. 

STEP  3)  a  =(4xl.049)/(2x3.1416x.3x3.75)  =  .5936 
STEP  4)  M=  206.37/1,160=  .18 

cc  =43.24  -  40.37  +  2.4  =  5.27  deg. 

Entering  into  curve  at  M= . 3  are  found  Cl=.87  and 
Cd  =  .015  . 

Y  =  arc  tan  (.015/. 87)  =  .988 


STEP  5)  assume  cti  =  5.85  deg. 

p2=  43.24  -  2x5.85  =  31.54  deg. 

j3^=  43.24  -  5.85  =  37.39  deg. 

STEP  6)  right  side  =((2xcos^37.39xcos.988(tan43.24- 

-tan31 . 54) )/ cos (37 . 39- . 988 )  = 

=  .51225 

left  side  =  .5936x.87  =.5164 

The  induced  angle  assumed  is  adequate  and  equal 


to  tti 

=  5 

.85  deg. 

STEP 

7) 

“cr  = 

5.27  -  5.85  = 

-.58  deg. 

STEP 

8) 

Clcr  = 

.41 

Cdcr  =• 

01 

Ter 

=  arc 

tan  (.01/. 

41)  =  1.397  deg. 

STEP 

9) 

Wm= 

150. 

34/COS37.39 

=  189.22  ft/sec 

Li  =  l/2x. 0023x189. 222x. 41x1. 049=17. 71  Ibf/ft 
STEP  10)  Ri  =  17.71  /  cosl.397  =  17.71  Ibf/ft 
Fi=  17.71  cos(37.39  -  1.397)  = 

=  14.33  Ibf/ft 
NEXT  STATION 
At  specified  station 
X=  .35 
c=  .958  ft 

|3  =  46.67  deg. 

STEP  1)  Cm  =  ( .2034(1- .35)+. 303)x337. 55  =  146.90  ft/sec 
U  =  2x3. 1416x. 35x3. 75x1200/60  =164.93  ft/sec 
W1  =  (146. 90^+164. 93^)^/^  =  220.86  ft/sec 
STEP  2)  Pi  =  arc  tan( 164 . 93 / 146 . 90 )  =  48.31  deg. 

STEP  3)  a  =(4x  .958)/(2x3.14l6x.35x3.75)  =  .4647 
STEP  4)  M=  220.86/1,160=  .19 

a=48.31  -  46.67  +  2.4  =  4.04  deg. 

Entering  into  curve  at  M= . 3  are  found  Cl=.83  and 
Cd  =  .015  . 

v=  arc  tan  (.015/. 83)  =1.033  deg. 


STEP  5)  assume  4.00  deg. 

^2-  48.31  -  2x4.00  =  40.31  deg. 

183=  48.31  -  4.00  =  44.31  deg. 

STEP  6)  right  side  =((2xcos^44.31xcosl.033(tan48.31- 

-tan40.31))/cos(44. 31-1. 033)  = 

=  .3864 

left  side  =  .4647x.83  =.3866 

The  induced  angle  assumed  is  adequate  and  equal 
to  «£=  4.00  deg. 

STEP  7)  arr-=  ^•0'^  -  ^-00  =  •0'^  deg. 


STEP  8)  Cl^^  =.455 


=.01 


=arc  tan  (.01/. 455)  =  1.259  deg. 

STEP  9)  Wm=  146.9/COS44.3  =  205.25  ft/sec 

Li  =  l/2x. 0023x205. 25^x.455x. 958=21. 11  Ibf/ft 
STEP  10)  Ri  =  21.11  /  cosl.259  =  21.11  Ibf/ft 
Fi=  21.11  cos(44.31  -  1.259)= 

=  15.43  Ibf/ft 

For  both  station  at  X  =  .3  and  X  =  .35 
STEP  11)  T=  (14. 33+15. 43)/2  x3.75x.05x  .325x3.75  = 

=  3.40  Ibf-ft 

After  that,  adding  all  station  torques  and  multiplying 
by  the  rotational  speed,  the  power  is  found. 

For  the  second  stage,  the  same  procedure  can  be 
followed  with  some  adjustments  as  mentioned  before  in 
chapter  3,  section  E, 


APPENDIX  C 

DESCRIPTION  OF  COMPUTER  PROGRAMS 


The  following  pages  contain  the  FORTRAN  PROGRAMS 
developed  to  evaluate  the  blade  performance  and  the  tunnel 
losses . 

The  LOSS  FORTRAN  is  a  program  to  evaluate  the  losses 
around  the  circuit  of  the  large  Academic  Wind  Tunnel  of  NPS . 
This  program  uses  one  subroutine  called  SKIN  to  calculate 
the  skin  friction  when  the  Reynolds  number  is  given;  it  does 
this  by  means  of  equation  2.1.  The  main  program  uses  the 
equations  for  Ko  and  Rey  for  each  section  of  the  tunnel, 
given  in  chapter  two.  The  meaning  of  variables  used  in  the 
program  are  given  at  beginning  of  the  main  program  and 
subroutine.  The  listing  of  the  program  is  presented  in 
Appendix  D. 

The  BORST  FORTRAN  is  a  program  used  to  evaluate  the 
performance  of  the  blades  of  both  stages.  This  program  uses 
the  procedure  presented  in  chapter  three  section  D,  to 
evaluate  the  torque  and  power  of  the  first  stage,  and 
section  E  to  evaluate  the  torque  and  power  of  the  second  ; 
it  uses  six  subroutines  as  follows:  FIRST  is  a  subroutine 
to  evaluate  the  parameters  of  the  blade  of  each  station  at  a 
specified  test  section  velocity,  for  the  first  stage. 
SECOND  is  a  subroutine  like  FIRST  but  it  evaluates  for  the 
second  stage.  CLNEW  is  a  subroutine  to  compute  the  lift 
coefficient  for  each  section  of  the  blade  using  Newton’s 
Forward  and  Backward  Interpolation  Formula,  with  the  data 
given  by  curves  of  Cl  vs  a  in  reference  12  .  CDNEW  is  a 
subroutine  to  compute  the  drag  coefficient  for  each  section 
using  the  same  formula  as  in  CLNEW,  with  data  given  by  the 
curve  of  Cd  vs  Cl  in  the  same  reference  as  in  CLNEW.  INDUC 
is  a  subroutine  used  to  compute  the  induced  angle  of  attack 


f  .if.it-. 


?S::a 


through  iterations  using  the  equation  3,4,  LIFT  is  a 
subroutine  to  compute  the  torque  and  power  for  second  stage 
of  blades  calling  all  the  subroutines  except  FIRST;  this 
subroutine  is  called  by  the  main  program  after  it  evaluates 
the  torque  and  power  of  the  first  stage. 

This  program  is  listed  in  Appendix  E  and  a  description 
of  parameters  appears  at  the  beginning  of  either  the  main 
program  or  the  subroutine. 
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COMPUTER  PROGRAM  '  LOSS 
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